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Foreword 


Highly pathogenic avian influenza (HPAI) caused by the H5N1 virus is a dangerous 
transboundary animal disease (TAD) that seriously impacts on poultry production, 
especially on smallholder farmers and commercial enterprises in Member States 
(MS) in Asia and Africa. With the added danger that this disease can become 
zoonotic, i.e. it has the potential to transfer to humans, it is essential that every effort 
is made to understand its epidemiology, particularly those factors that determine the 
dissemination of the virus across countries and continents. One area in which there 
is an urgent need for clarity is the role that migratory wild birds might play in the 
spread of HPAI to domestic poultry. This is a matter of international importance 
since millions of wild birds fly from areas where HPAI is endemic to all corners of 
the globe, potentially spreading the disease wherever they land. 

MS are acutely aware of the need to develop technologies that allow disease 
detection before clinical signs of infection became apparent. Under the aegis of the 
project “Reducing risk from transboundary animal diseases (TADs) and those of 
zoonotic importance”, the IAEA, through the Animal Production and Health Section 
of the Joint FAO/IAEA Division, has supported MS in acquiring sensitive, specific 
and rapid nuclear, nuclear-related and nuclear-associated diagnostic tests for 
HPAI. Now that MS laboratories have the required technologies to develop appro¬ 
priate surveillance strategies, it will be possible to enlarge the scope of studies on 
HPAI to include verification of the role of the wild bird hosts. Consequently, the 
Animal Production and Health Section (APHS) has initiated a new Coordinated 
Research Project (CRP) that uses innovative technologies to sample wild birds and 
their habitats to detect virus and also introduce stable isotope analysis (SIA) of bird 
tissues to determine their place of origin. 

Migration patterns of wild birds have been studied for many years using a variety 
of techniques, including recent various electronic tracking methods. However, more 
effective and cheaper options will be required if specific information is needed on a 
particular bird species, especially infected wild birds. The envisaged approach is to 
determine migrant origins by determining the stable isotope “signature” of the bird 
and match it to the most likely environmental habitat. Stable isotopes pose no health 
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or environmental risks since they are not radioactive and are present throughout the 
environment. However only a small number of them, namely, H, O, C, N and S, will 
be useful for tracking the origin of wild birds. The technologies used to infer migra¬ 
tory connectivity in migrant birds and for detection have become increasingly 
sophisticated, as have the models for identifying habitats. It is now felt that such 
technologies are mature enough to inform scientists engaged in studies of HPAI on 
so far unanswered questions concerning the biology and epidemiology of HPAI in 
wild birds. 

The purpose of this publication is to make scientists, particularly those involved 
in IAEA CRP and TC projects, aware of the potential of SIA by providing an intro¬ 
duction of their application and technical requirements needed to facilitate studies. 

Animal Production and Health Section Gerrit Viljoen 

Joint FAO/IAEA Division for Nuclear Applications 

in Food and Agriculture, Department of Nuclear 

Sciences and Applications 

International Atomic Energy Agency 

Vienna International Centre 

Vienna, Austria 
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Chapter 1 

General Introduction 


1.1 Background Analysis 

Birds are amongst the few terrestrial vertebrates that share with humans the pecu¬ 
liarity of travelling in a few hours across national and intercontinental borders 
(Hobson 2002). The record for distance covered in a single year belongs to the 
arctic tern, which travels over 50,000 km between Antarctica and northern 
Scandinavia. Overall, billions of birds travel between continents twice a year in 
only a few weeks. Migration is critical in the life cycle of a bird, and without this 
annual journey many birds would not be able to raise their young. More than 5000 
species of birds manage annual round-trip migrations of thousands of miles, often 
returning to the exact same nesting and wintering locations from year to year. Birds 
migrate to find the richest, most abundant food sources that will provide adequate 
energy to nurture young birds. If no birds migrated, competition for adequate food 
during breeding seasons would be fierce and many birds would starve. Instead, 
birds have evolved different migration patterns, times and routes to give themselves 
and their offspring the greatest chance of survival. Birds gauge the changing of the 
seasons based on light level from the angle of the sun in the sky and the amount of 
daily light. When the timing is right for their migrating needs, they will begin their 
journey. Several minor factors can affect the precise day any bird species begins its 
migration, including available food supplies, poor weather or storms and air tem¬ 
peratures and wind patterns. While these factors may affect migration by a day or 
two, most bird species follow precise migration calendars. Those calendars vary 
widely for different species, however, and while autumn and spring are peak migra¬ 
tion periods when many birds are on the move, migration is actually an ongoing 
process and at any time of the year, there are always birds at some stage of their 
journeys. The distance the birds must fly, the length of time it takes to mate and 
produce a healthy brood, the amount of parental nurturing young birds receive and 
the location of birds’ breeding and wintering grounds all affect when any one spe¬ 
cies migrates to stay alive. 
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During these yearly migrations, birds have the potential of dispersing microor¬ 
ganisms that can be dangerous for human as well as animal health. For instance, 
birds are believed to be responsible for the wide geographic distribution of various 
pathogens, including viruses [e.g., West Nile, Sindbis, Highly Pathogenic Avian 
Influenza (HPAI), Newcastle disease], bacteria (e.g., Borrelia, Mycobacterium, 
Salmonella ), and protozoa (e.g., Cryptosporidium). An insight into the ecology of 
bird populations is necessary to understand fully the epidemiology of bird-associated 
emerging diseases. Furthermore, data about avian movements might be used to 
improve disease surveillance or to adapt preventive measures. However, the links 
between bird ecology and livestock and human disease have yet to be completely 
understood so there is a need to increase knowledge of avian migration patterns and 
infectious diseases to help predict future outbreaks of emerging diseases. 

Animal influenza viruses threaten animal health, livestock productivity and food 
security in poor countries, but they can also evolve into dangerous human patho¬ 
gens. This has been seen with the emergence of HPAI. Its main impact has been on 
domesticated poultry, with over 300 million birds killed or destroyed, but of consid¬ 
erable public health concern is the transmission of the virus from birds to humans, 
resulting in over 500 occurrences of disease in which over half of the infected indi¬ 
viduals have died. 

The threat from animal influenza viruses makes it essential for animal health 
professionals to take the lead in detecting and monitoring the occurrence of the 
viruses and sharing the information with the international community. Since avian 
influenza appears to be associated with migratory bird movements (over 1000 
reported AI outbreaks since 2006, involving 25 species of wild birds in EU alone) 
surveillance would need to focus on detection of HPAI in both wild birds as well as 
domesticated poultry and it will also be necessary to establish the migratory path¬ 
ways of wild birds to increase the capacity to assess their risk in spreading the virus. 

In March 2010 the International Scientific Task Force on Avian Influenza and 
Wild Birds, led by FAO and the United Nations Environment Programme - 
Convention on Migratory Species, reported that waning attention to HPAI was 
reducing opportunities for surveillance and research, thereby affecting efforts to 
understand the epidemiology of the disease. The disease continues to be a major 
problem in Egypt and parts of Asia and outbreaks have occurred in poultry in 
Romania and in wild birds in Russia, China and Mongolia (Fig. 1.1). One of the 
most important issues that need addressing is the surveillance of wild bird popula¬ 
tions to improve understanding of the role that they play in the dissemination of 
infection. Although various methods have been used to track the migration of birds 
there is increasing interest in utilizing methodologies that would enable tracing of 
migratory movement based on the birds’ stable isotope signatures. A small number 
of stable isotopes are involved in important biological and ecological processes and 
there is a strong correlation between levels of these isotopes in the environment and 
the concentration of the same isotopes in avian tissues. Of most interest are stable 
hydrogen and oxygen ratios in tissues that accurately reflect those in lakes, rivers, 
oceans and in groundwater, along the wild bird fly ways. Using stable isotopes to 
characterize a population involves examining the isotopic signatures of a few 
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Fig. 1.1 Characteristics of the reported outbreaks of AI in domesticated poultry and wild migra¬ 
tory birds in the period from 1 January 2011 until 13 May 2015 (Source: FAO 2015). 

Map A symbols: red circle [#]=H5 AI; blue square |D| = H7 AI; green triangle [A] = H9 AI; 
grey star [lV]; 

Map B symbols: AI outbreaks reported during 2011 ( grey star [h]), 2012 (, green square [□]), 
2013 {blue square [^]), 2014 {red triangle [A]) and 2015 {yellow circle [O]); 

Map C symbols: AI outbreaks in wild {red circle [#]) and captive birds {blue triangle [A]) 


individuals that are representative of the entire population. The hydrogen and oxy¬ 
gen isotope composition of environmental water varies spatially across the globe 
and because it is a constituent of many biosynthetic pathways, the isotopes’ pres- 
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ence is relayed to animal tissues, providing the means to link data on groundwater 
isoscapes with isotope levels in biological tissues such as feathers. This isotope data 
would reveal migration patterns and enable identification of the breeding areas of 
birds sampled at non-breeding grounds and disease outbreak sites. 


1.2 Migratory Birds and HPAI 

To date, only a small number of migratory birds have been tracked by satellite trans¬ 
mitters to establish links with disease outbreaks in domestic poultry, with little evi¬ 
dence so far of direct correlation with HPAI. Furthermore, surveillance of 750,000 
“healthy” wild birds has not revealed many infected individuals. Circumstantial 
evidence suggesting that the spread of the H5N1 virus to new areas can be facilitated 
by migratory wild birds has come from studies on ducks in Asia marked with satel¬ 
lite transmitters (Yamaguchi et al. 2010) that were tracked during an outbreak of 
highly pathogenic H5N1 avian influenza virus. 

The satellite transmitters were attached to northern pintail ducks several months 
before the H5N1 virus was discovered in dead and dying whooper swans at wet¬ 
lands in Japan. Twelve percent of marked pintails used the same wetlands as 
infected swans and the pintail ducks were present at those sites on dates the virus 
was discovered in swans. During the first week after they become infected with 
H5N1 virus, ducks such as pintails can shed the virus orally or in their faeces, con¬ 
tributing to the virus’ spread. Some of the marked pintails migrated 700 miles 
within 4 days of leaving the outbreak sites; marked pintails ultimately migrated 
more than 2000 miles to nesting areas in eastern Russia. The discovery that northern 
pintails made long-distance migrations during the period when an infected duck 
would likely shed the virus offers insight into how H5N1 could be spread by wild 
birds across large areas (Yamaguchi et al. 2010). 

In another study, the movements of bar-headed geese marked with GPS satellite 
transmitters at Qinghai Lake, China were traced in relation to virus outbreaks and 
disease risk factors. A previously undocumented migratory pathway between 
Qinghai Lake and the Lhasa Valley of Tibet where 93 % of the 29 marked geese 
overwintered was discovered. From 2003 to 2009, 16 outbreaks in poultry or wild 
birds were confirmed on the Qinghai-Tibet Plateau, and the majority were located 
within the migratory pathway of the geese. Spatial and temporal concordance 
between goose movements and three potential H5N1 virus sources (poultry farms, 
a captive bar-headed goose facility, and H5N1 outbreak locations) indicated that 
ample opportunities existed for virus spill over and infection of migratory geese on 
the wintering grounds. Their potential as vectors of H5N1 was supported by rapid 
migration movements of some geese and genetic relatedness of H5N1 virus isolated 
from geese in Tibet and Qinghai Lake. This study was the first to compare phyloge¬ 
netics of the virus with spatial ecology of its host, and the combined results suggest 
that the geese play a role in the spread of H5N1 in this region (Prosser et al. 2006; 
Zhou et al. 2011). 
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Defining the migratory behaviour of animals is fundamental to understanding 
their evolution and life history; much of the effort to study migration in birds, has 
been driven by the needs for conservation in changing habitats. The advances in 
understanding of migration, particularly in North America and Europe have been 
aided considerably by stable isotope techniques. The methods have not yet been 
applied to study the migratory birds and the dissemination of HPAI, but this is likely 
to be remedied in the near future. This manual provides background information on 
Stable Isotope Analysis (SIA) for scientists investigating HPAI from Member States 
engaged in Collaborative Research Projects or Technical Contract Projects and wish 
to understand how the disease is disseminated to provide information that would 
enable better risk assessment and more effective control and prevention measures. 

In order to follow the global distribution of animal diseases, as well as to enable 
rapid information sharing, FAO has established an Animal Disease Information 
System (EMPRES-i), (FAO 2015) which enables for filtering disease records 
according to different predefined criteria. Specifically for AI, record-sets can be 
generated, which enable for correlation according to the species of origin of the 
outbreak (domesticated or wild birds), period of the observation/reporting of the 
outbreak and the location of the outbreak. The data on the evolution of the AI glob¬ 
ally, obtained from EMPRES-i system over the period of 5 years (January 2011- 
May 2015) is shown on Fig. 1.1 and Table 1.1. 


1.3 Using SIA to Understand the Dissemination 
of HPAI - The Way Ahead! 

Most elements consist of one or more stable isotopes - elements having the same 
number of protons, but differing in the numbers of neutrons. Stable isotopes are 
those isotopes of an element that do not decay through radioactive processes over 
time. For instance, the element carbon (C) exists as two stable isotopes, 12 C and 13 C, 
and the element hydrogen (H) exist as two stable isotopes, ! H and 2 H. 

Stable isotope contents are expressed in ‘delta’ notation as 5 values in parts per 
thousand (%c), where: 


8R%c = (Rs / RStd -1) x 1000 

and Rs and RStd are the ratios of the heavy to light isotope (e.g. 13 C/ 12 C) in the 
sample and the standard, respectively. The stable isotope ratios of hydrogen, carbon, 
nitrogen, oxygen and sulphur are denoted in delta notation as 5 2 H, 5 13 C, 5 15 N, 5 18 0 
and 5 34 S, respectively. R values have been carefully measured for internationally 
recognized standards. The standard used for both H and O is Standard Mean Ocean 
Water (SMOW), where (WH) standard is 0.0001558 and 18 0/ 16 0 is 0.0020052. 
The original SMOW standard is no longer available and has been replaced by a new 
International Atomic Energy Agency (IAEA) standard, VSMOW. The international 
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carbon standard is the PDB, where ( 13 C/ 12 C) is 0.0112372 and is based on a belem- 
nite from the Pee Dee Formation. As with SMOW, the original PDB standard is no 
longer available, but the IAEA provides Vienna-PDB with a similar R value. 
Atmospheric nitrogen is the internationally recognized standard with an R value 
( 15 N/ 14 N) of 0.0036765. Lastly, the internationally recognized standard for sulfur is 
CDT, the Canyon Diablo Troilite, with a value ( 34 S/ 32 S) of 0.0450045. Typically, 
during most stable isotope analyses, investigators would not use IAEA standards on 
a routine basis. Instead, laboratories establish secondary reference materials to use 
each day that are traceable to IAEA standards and that bracket the range of isotope 
ratio values anticipated for the samples. 

Although ecologists refer to “isotope signatures or isotope fingerprints” the val¬ 
ues obtained for a sample do not provide a unique fingerprint but a distinctive pro¬ 
file; a more neutral notation is isotope value, rather than signatures. 

The 5 notation is derived as follows: 


SH(% 0 ) = 


r Isotopic ratio sample 
v Isotopic ratio standard 




-1 


xlOOO 


J 


The right side of the equation is the measure of the light to heavy isotope ( 2 H/ ! H) 
and that ratio is multiplied by 1000 transform the values into whole numbers. The 
isotope reference points were established many years ago and results reported can 
be negative or positive (%c) relative to the accepted international standard. Thus, a 
5H value of +150 (%c) means the sample has 150 parts per thousand (15 %) more 
deuterium in it than the standard; while if the value were negative it would be 150 
times less deuterium than the standard. Primary reference material is limited in 
quantity so laboratories tend to use local standards calibrated against a reference 
standard. The primary isotopic reference standard for hydrogen and oxygen is 
Vienna Standard Mean Ocean Water (VSMOW). 

It will be necessary to develop suitable procedures for monitoring wild birds to 
detect HPAI, determine their origin and migratory pathway to estimate the duration 
of their stay at the place of capture. This will require the application of detailed 
standard operating procedures for feather and tissue collection, data recording, sam¬ 
ple storage, and designation of a laboratory to carry out the SIA. Since tissue sam¬ 
ples could come from birds suspected of carrying avian influenza virus it will be 
necessary to install safeguards to prevent possible infection. Soft tissues could be 
freeze-dried and both these and keratinaceous samples heated to 100 °C for 20 min 
to destroy any viruses. Alternative methods would be to irradiate specimens at a 
central location before making them available for study. Various innovative meth¬ 
ods for virus detection should also be considered, including environmental sam¬ 
pling to detect H5N1 virus in water, or faeces (Khalenkov et al. 2008; Cheung et al. 
2009; Dovas et al. 2010) in locations where wild birds congregate. Other proce¬ 
dures that will be useful could include identification of bird species using DNA 
barcoding from feather or faeces (Lee et al. 2010) to provide better links to identify¬ 
ing those birds that are truly implicated in dissemination of HPAI. 
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There are already available data on variation in stable isotope ratios across the 
globe and how this is reflected in tissue samples from birds inhabiting different 
regions. It might be necessary however to collect environmental samples from dif¬ 
ferent locations if there is not sufficient data already available. Collection of wild 
birds is a specialized procedure and it will be necessary to liaise with wildlife groups 
skilled in this task, who are also able to provide basic data on the ecology of differ¬ 
ent bird species. The FAO Wildlife Unit can provide inputs in this area through its 
links with various wildlife organizations in order to carry out sampling with greater 
efficiency. 

Consideration needs to be given to which isotopes to analyse. While stable 
hydrogen can be used to obtain information on the wider geographical location, 
further details might be inferred on the local habitat by analysing carbon, nitrogen 
and sulphur. Samples from keratinaceous tissues e.g. feathers and claws, will pro¬ 
vide an isotope ratio of the place when they were grown. In blood, stable isotope 
ratios could be used to determine timing of arrival on breeding or wintering grounds. 
Isotope ratios in the blood relate to those in the current environment hence a differ¬ 
ence between isotope ratios in blood and the environment would indicate newly 
arrived birds compared with birds that had been staying in a particular location for 
some time. If infected birds were to be found it should be possible to establish 
where the disease was acquired - at the site or in a previous location. By under¬ 
standing migratory movements it will be possible to predict risk and derive models 
to show the spread of HPAI and identify areas that pose a significant potential for 
being locations where the dissemination of HPAI is greater due to the particular 
congregation of wild birds, H5N1 virus and domestic poultry. 

Open Access This chapter is distributed under the terms of the Creative Commons Attribution- 
Noncommercial 2.5 License (http://creativecommons.Org/licenses/by-nc/2.5/) which permits any 
noncommercial use, distribution, and reproduction in any medium, provided the original author(s) 
and source are credited. 

The images or other third party material in this chapter are included in the work’s Creative 
Commons license, unless indicated otherwise in the credit line; if such material is not included 
in the work’s Creative Commons license and the respective action is not permitted by statutory 
regulation, users will need to obtain permission from the license holder to duplicate, adapt or 
reproduce the material. 


Chapter 2 

Animal Migration Tracking Methods 


2.1 Extrinsic Markers 

The migration of birds has been studied for many years relying mostly on extrinsic 
passive markers attached to individual animals at the point of capture, with the 
expectation that a proportion of the marked individuals will then be identified in 
another location at a different point in time. Over the past 100 years, the most wide¬ 
spread approach has been through the application of markers such as leg bands, 
neck collars, or dyes. Many millions of birds have been tagged in this way but 
although this method has provided insights into migration, for the vast majority of 
bird species examined the recovery rate is low. An alternative to these simple devices 
is to use miniature transmitting devices - radio transmitters, radar and satellite 
tracking - that serve as active markers and are small enough (<0.5 g) to be attached 
to even small birds or mammals. The location of the marked animal can be inferred 
by tracing the individual using a receiver, or by triangulation using several receivers. 
Since the devices are miniaturized their range and battery life are restricted and they 
can provide information over only a few kilometers. Radar technology has also 
made useful contributions to studies on migration since it can provide information 
on animal movements over considerable distances, but since radar installations are 
fixed it is not possible to trace movements over the whole spectrum of migration 
routes used by birds. The most significant advances in tracking migratory animals 
have come from the use of satellite transmitters that allow highly accurate position¬ 
ing of individual animals (Hiroyoshi and Pierre 2005; Whitworth et al. 2007; http:// 
www.fao.org/avianflu/en/wildlife/sat_telemetry.htm). Much of the globe is covered 
by satellites so that animals can be monitored over thousands of kilometers. The 
technique can only be used on relatively large animals as the weight of the smallest 
transmitters is approximately 10 g, restricting their use to an animal weighing about 
250 g, thereby excluding 80 % of the world’s birds and 70 % of the mammals. With 
the exception of satellite transmitters all extrinsic markers require that individuals 
be recaptured, re-sighted or move within a detector’s range at some time after initial 
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marking. The probability of recapture depends on the number of observers, the 
regions and habitats and the chances of success are low. In addition, extrinsic meth¬ 
ods tend to be biased towards regions with a high likelihood of mark-recapture 
(Hobson et al. 2004). A fundamental flaw in the use of an extrinsic marker is that it 
provides information only on the marked individuals. Geolocators and satellite 
tracking rely on small sample sizes and the devices may affect the behavior of the 
marked bird (Stutchbury et al. 2009). Extrapolating the findings from one individual 
to a whole population depends on how representative the marked individuals are. A 
single recovery or satellite track may not reveal what the population is doing. 


2.2 Intrinsic Markers 

The major advantage of an intrinsic marker is that it is not necessary to capture and 
mark an individual in advance, and every capture provides information on the origin 
of that animal. There is therefore, less bias than with extrinsic markers, where the 
location of the origin is dependent on accessibility; for example, remote areas might 
not be easily included in any migration study. Man-made contaminants such as 
dioxins and methyl mercury or other heavy metals could provide a means of infer¬ 
ring migratory origins as exposure to them varies throughout the world. Exposure 
of migratory animals to parasites or pathogens also varies geographically and could 
be used for determining movements but little research has been conducted as to how 
such markers might assist in deciphering migration patterns. In a study of Sharp- 
shinned Hawks ( Accipiter striatus) in New Mexico, stable-hydrogen isotope ratios 
of feathers were estimated and blood was collected to quantify the prevalence and 
intensity of haematozoan infection in the birds (Smith et al. 2004). Twenty four 
percent of the birds were infected with Leucocytozoon toddi , 37 % with 
Haemoproteus elani , and 5 % with H. janovyi. This was the first documented occur¬ 
rence of H. janovyi in North America in these hawks and the stable-hydrogen iso¬ 
tope analyses indicated that only birds originating from south-western North 
America harboured H. janovyi which could be of significance in regard to biogeog¬ 
raphy of that parasite species. The isotope studies also showed that birds infected 
with Haemoproteus spp were more widespread geographically than those with 
Leucocytozoon species. Another type of intrinsic markers that could potentially 
track animal movement geographically is the trace elements. Animals acquire quite 
distinctive chemical profiles relative to a particular geographical location and retain 
that profile when they move to another area. Trace element profiles have been used 
to determine if populations of breeding birds originate from different winter loca¬ 
tions (Sze'p et al. 2008). However, there is little information on how trace elements 
vary over the globe, so it is of only limited potential. Stable isotope tools to trace 
animal migration have been used extensively in ecological studies but have not yet 
been utilized in work on the dissemination of transboundary animal diseases. There 
is potential for greater understanding the epidemiology of HPAI in relation to the 
migration of wild waterfowl by using stable isotope studies to identify the origins and 
stopover points of birds suspected of carrying HPAI. A summary of the advantages 
and disadvantages of extrinsic and intrinsic methods is given in Tables 1.1 and 2.1. 
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Table 2.1 Extrinsic markers for tracking animal migration 


Technique 

Advantages 

Disadvantages 

Extrinsic 

1. Can apply to a 
broad range of 
animals 

1. Requires capture/recapture 

2. High spatial 
resolution 

2. Biased towards initial 
captured population 

Phenotypic variation 

1. Inexpensive 

1. Low resolution 

2. Can be used for 
historical specimen 

2. Not applicable to all 
species 

Banding/Marking 

1. Inexpensive 

1. Low recovery rates <0.5 % 

2. Gives 

information on start 
and end of 
migration 

2. Data must be acquired 
over many years 

3. Only few banding stations 
across the world 

4. Marking and recovery are 
biased towards the banding 
station location 

Radio-transmitter 

1. Gives precise 
locations 

1. Low range 

2. Gives precise 
trajectory 

2. Expensive 

3. Transmitters might affect 
behaviour 

Satellite transmitters 

1. Precise animal 
trajectory 

1. Expensive 

2. Only for use on large 
animals 

3. Transmitters might affect 
behaviour 

International Cooperation for Animal 
Research Using Space ICARUS (www. 
Icaruslnitiative.org) 

1. Transmitter 
allows many 
different species to 
be tracked 

1. High start-up investment 

2. Can track 
individuals over the 
globe 

2. Technology not proven 

Satellite tracking to determine the 
spread the spread of infectious diseases 

3. Inexpensive after 
start-up costs 

3. Transmitters might affect 
behaviour 

Passive radar 

1. Coverage over 
large geographical 
region 

1. Coverage only from 
existing stations or portable 
instruments 

2. Inexpensive 

2. Poor ability to determine 
species and individual 
movements 

3. Individuals need 
not be captured 

Transponders 

1. Small size 

1. Requires external 
activation 

2. Low range 

3. Coverage only form 
existing stations or portable 
instruments 


(continued) 
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Table 2.1 (continued) 


Technique 

Advantages 

Disadvantages 

Geolocation tags 

1. Produces animal 
trajectory 

1. Individuals must be 
captured to download data 

2. Low weight 

2. Accuracy relative to 
satellite tags low 

3. Inexpensive 
relative to satellite 
tags 


Studies of avian migration and population structure take advantage of geographi¬ 
cal variation in stable isotope ratios, in particular the ratio of the two stable isotopes 
of hydrogen in precipitation. In North America, for instance, the ratio of growing 
season precipitation decreases from -20%c VSMOW (Vienna Standard Mean 
Ocean Water) in Florida to -140%c VSMOW in north-western Canada (McKechnie 
2004). Because the deuterium of precipitation is incorporated into avian tissues via 
food webs and feathers are metabolically inert once grown, the deuterium values of 
feathers reflect the latitude at which they were grown. In species that moult prior to 
migration, feather deuterium can be used to identify the latitudinal origin of indi¬ 
viduals, and hence link breeding and wintering grounds. Isotopic data can also be 
combined with information on capture dates or ring recovery. Stable isotope studies 
require a detailed knowledge of a species’ moult schedule, and the assumption that 
the isotope signature of feathers reflects that of food webs in the region where the 
feather is grown must be verified. The technique has the best potential where there 
is a good spatial variation in the environmental locations in which the particular bird 
species is found. 

The basic requirements in the use of stable isotopes are that: 

• The characteristics of the environment through which the animal of interest 
moves are known. The term used to describe the mapping of large-scale and 
spatiotemporal distributions of stable isotope ratios in the environment and ani¬ 
mals is an isoscape or isotopic landscape (Bowen and West 2008). 

• Isotopic values in animal tissues can be discriminated from baseline isoscape 
values, 

• The time period for acquisition into a particular animal tissue is understood. For 
instance, turnover of isotope will be greater in metabolically active tissues than 
in tissues with a slow, or no, turnover e.g. feathers (Table 2.2). 

Just a handful of these stable isotopes are of any practical interest to studies on 
animal migration. These are the “light isotopes” of H, O, C, N, and S that are pres¬ 
ent in all components of the biosphere (plants and animals), the hydrosphere (water) 
and the atmosphere (gaseous 0 2 , N 2 and H 2 0). These five elements comprise the 
bulk of all animal tissues (Table 2.3). These elements and their isotopes circulate in 
the biosphere to produce characteristic isotope distributions globally. Large pools of 
these elements provide stability in the overall isotope circulation e.g. in the ocean, 
for hydrogen and oxygen isotopes, inorganic carbon pool in the ocean for carbon 
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Table 2.2 Intrinsic markers for tracking animal migration stable isotopes 


Technique 

Advantages 

Disadvantages 

Intrinsic 

1. Not biased to initial capture 
population 

1. Biased to final capture population 

2. Less labour intensive than 
most extrinsic methods 

2. Lower resolution than extrinsic 
methods 

Contaminants 

1. Potentially high spatial 
resolution 

1. Lack of distribution maps 

2. Transport of contaminants could 
give unreliable geographical signal 

Parasites/Genetics 

1. Several possible markers 

1. Species specific 

2. Low resolution 

Trace elements 

1. Measurement of a large 
number of elements 

1. Lack of distribution maps 

2. Potentially high spatial 
resolution 

2. Expensive 

3. Requires more sample tissue 

4. Requires tissue that is 
metabolically inactive after growth 

5. Spatial resolution could be too 
high 

6. Some elements may be integrated 
into inactive tissue after growth is 
complete 

Stable isotopes 

1. Inexpensive 

1. Low resolution 

2. Not species- or 
taxon-specific 

2. No base maps 

3. Several isotopes can be 
combined to increase spatial 
resolution 

3. Ideal tissue is metabolically 
inactive after growth 

4. Turnover rate of elements in 
active tissue is unknown 

5. Interpretation may be complicated 
by animals’ physiology 


isotopes, and sulphate in the sea for sulphur isotopes and the atmospheric reservoir 
for nitrogen. The presence of the stable isotopes varies widely in nature. All of the 
light isotopes have a common or abundant form, e.g. Tf; 99.985 % and a more rare 
“heavier form”, 2 H; 0.015 % (Table 2.3). The abundance ratios of these isotopes 
vary because of physical and chemical processes they undergo in nature and it is 
these variations that enable the use of stable isotopes in tracing migration. As it is 
difficult to measure the precise concentrations of the isotopes in samples, measuring 
the relative differences in isotopic ratios between a sample and a reference by means of 
a mass spectrometer is the way in which data are acquired. Since gas source isotope 
ratio mass spectrometers (IRMS) are used to measure light isotopes it is not possible 
to measure isotope ratios directly from a tissue sample - feather, blood, muscle, 
claw, hair - instead, the sample must first be combusted to form a gaseous analyte 
in an elemental analyzer, a gas chromatograph or a laser and the resulting gas can 
then be used to measure the isotopic ratios relative to a calibrated reference gas of 
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Table 2.3 Dry weight % abundance of light stable isotope ratios of interest in determining 
migratory connectivity in tissues 


Element 

Weight (%) 

Isotope ratios 

5 Range (0/00) 

Mass required (mg) 

Carbon 

30-40 

13 C / 12 C 

05 to -65 

0.2-1.5 

Oxygen 

27-40 

i8 o/ 19 o 

+10 to +30 

0.2-00.5 

Nitrogen 

12-19 

15 N / 14 N 

—2 to +25 

0.5-1.5 

Hydrogen 

6-8 

2 H/‘H 

-250 to +90 

0.1-0.4 

Sulphur 

5-20 

34 S/ 32 S 

-20 to +30 

1-2 


Table 2.4 Average terrestrial 
abundances of the stable 
isotopes of major elements of 
interest in ecological studies 
the carbon cycle 


Element 

Isotope 

Abundance (%) 

Hydrogen 

*H 

99.985 

2 H 

0.015 

Carbon 

12 C 

98.89 

13 c 

1.11 

Nitrogen 

14 N 

99.63 

15 N 

0.37 

Oxygen 

i6 o 

99.759 

17 o 

0.037 

i8 o 

0.204 

Sulphur 

32 S 

95.00 

33 S 

0.76 

34 s 

4.22 

36 s 

0.014 


the same type. Development of IRMS equipment has led to automatic sample pro¬ 
cessing to obtain multiple isotope assays from a single sample, on decreasing sam¬ 
ple size and improving throughput rates. IRMS can measure isotopic ratio differences 
to the sixth decimal place or ±0.01 % (Table 2.4). 

The carbon cycle involves active exchanges of C0 2 among the atmosphere, ter¬ 
restrial ecosystems and the surface ocean. The 13 C value of atmospheric C0 2 is 
decreasing in response to inputs of 13 C depleted C0 2 from fossil fuel plus biomass 
burning and decomposition. Over the past 100 years the decrease may have been 
almost l%o, from about -l%c to -8 %c. 

Depending on how plants’ photosynthesis process is materialized, they are clas¬ 
sified in two large groups, C3 and C4, with very different values for 5 13 C. In the C3 
group, the first photo synthesized organic compound has 3 atoms of carbon while in 
group C4, there are 4. Most plants (85 %) (E.g. trees and crops) follow the C3 pho¬ 
tosynthesis pathway and have lower values of 5 13 C, between -22%c and -30%c. The 
remaining 15 % of the plants are of type C4. The majority are tropical herbs and 
have high values of 5 13 C, between -10%c and -14%c (http://homepage.mac.com/ 
uriarte/carbon 13 .html) . 

Carbon uptake by the dominant C3 plants on land involves a net fractionation of 
about 20%c between the atmospheric C0 2 and plant biomass (-28%c). Carbon 
uptake by C4 plants, mainly tropical and salt grasses, involves a small net fraction- 
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ation of about 5 %c. Soil organic matter globally contains several-fold more carbon 
than either the atmosphere or living plant biomass and is similar or slightly enriched 
in 13 C in comparison with the dominant vegetation. The exchange of C0 2 between 
the atmosphere and the surface of the ocean involves an equilibrium chemical frac¬ 
tionation between atmospheric C0 2 (-8 %o) and the total C0 2 in surface ocean water. 


2.2.1 The Nitrogen Cycle 

Most nitrogen in the biosphere is present as N 2 gas in the atmosphere. This massive 
reservoir is well mixed with an isotope composition that is essentially constant at 
Woo. Nitrogen in most other parts of the biosphere also has an isotope composition 
near the 0 %c value, from -10 to +10%c, primarily because the rate of nitrogen sup¬ 
ply often limits reactions such as plant growth and bacterial mineralization. Under 
these conditions all available nitrogen can be consumed, without regard to isotope 
content and with no overall isotope fractionation. Thus, slow rates of N supply and 
limiting amounts of substrate N are often important for understanding nitrogen iso¬ 
tope distributions. Some cumulative and large fractionations do occur in the nitro¬ 
gen cycle. Lakes appear more variable in isotope composition than the large world 
ocean. Large isotope contrasts might be expected between lakes in which primary 
production is limited by N (little fractionation by phytoplankton) versus P (abun¬ 
dant large possible fractionations during N uptake by phytoplankton). Where 
phytoplankton have different 15 N values than terrestrial vegetation, the nitrogen iso¬ 
topes may function as source markers for autochthonous and allochthonous organic 
matter. 


2.2.2 The Sulphur Cycle 

Sulphate in the ocean is a large well-mixed sulphur reservoir whose isotope compo¬ 
sition is 21%c heavier than primordial sulphur in the earth and solar system at large. 
Fixation of sulphate by phytoplankton occurs with a small isotope effect, but sul¬ 
phate reduction in marine sediment occurs with a large effect of 30-70%c. Over 
geological time, and partially in response to global-scale fluctuations in sulphate 
reduction activities, the 34 S values of oceanic sulphate have varied from about +10 
to +33%c. Uplift and preservation of marine sedimentary sulphides and sulphate- 
containing evaporates on land have produced a patchwork of sulphur in terrestrial 
environments, each with different 34 S values for bedrock sulphur. Thus, large 34 S 
ranges must be assigned in general sulphur cycle diagrams. In spite of this, conti¬ 
nental vegetation seems to average near +2 to +6 %c over large areas and is quite 
distinct from the +~17 to +21%c values of marine plankton and seaweeds. The sta¬ 
ble isotope composition of sulphur entering the atmosphere can also be quite 
variable. 
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2.2.3 The Oxygen Cycle 

There are three oxygen isotopes that act as tracers when the many common oxygen- 
containing molecules circulate in the biosphere. The water cycle controls much of 
the oxygen dynamics and oxygen isotope dynamics. Evaporation and condensation 
result in predictable variations in isotope compositions of water that are now rou¬ 
tinely tracked at regional and global levels (http://isohis.iaea.org/userupdate/water- 
loo/index.html, http://www.waterisotopes.org/ and http://ecophys.biology.utah. 
edu/labfolks/gbowen/pages/Isomaps.html#IAEA). 

Oxygen isotope studies with animals have focused on determining which local 
sources of water are used. The degree to which food influences 18 0 variations in 
animals has not been determined fully. 


2.2.4 The Hydrogen Cycle 

Much of the hydrogen cycle involves water, with various processes in the water 
cycle leading to characteristic, large-scale geographic patterns of hydrogen isotopes 
in water. Ocean water is the main reservoir of hydrogen in the biosphere and the 
standard reference material “standard mean ocean water” for hydrogen isotope 
measurement. The isotope composition of ocean water is a good starting point for 
following isotope dynamics in the hydrological cycle. The transitions between liq¬ 
uid water and water vapour during evaporation and condensation involve kinetic 
and equilibrium reactions with isotope fractionation. Water vapour evaporating 
from the sea has 2 H values of -10 to -20%c, and as this process reverses during 
condensation and formation of rain and snow, this trend towards lower atmospheric 
2 H values is amplified. As water vapour moves inland and up mountains, it progres¬ 
sively loses moisture and 2 H values decline further. These processes can be ampli¬ 
fied yet again in colder regions where low temperatures promote stronger 
fractionations between vapour and condensate. A combination of high elevation 
and low temperature can result in 2 H values of -200 to -400%c for water in high- 
elevation glaciers and for snow in Polar Regions. In less dramatic examples, large 
rivers fed by snowmelt that have continental origins can have much lower 2 H values 
than coastal marine waters. This makes 2 H source signals valuable tracers in coastal 
estuaries and floodplains linked to these rivers. Isotope hydrology studies often con¬ 
sider the water isotopes (hydrogen and oxygen isotopes) as markers for water 
sources and water circulation. Global-scale maps and animations of hydrogen and 
oxygen isotope variations in water are available on the Web (http://isohis.iaea.org/ 
userupdate/waterloo/index.html or www.waterisotopes.org). 

Analytical advances are making it easier to investigate the origins and cycling of 
hydrogen bound in organic matter. About 10-20 % of hydrogen in organic materials 
is exchangeable with water vapour present in normal laboratory air, but this 
exchange effect is understood and can be corrected for during routine analysis. 
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Hydrogen in animal tissues can be divided into three main pools: hydrogen derived 
from dietary sources, hydrogen from drinking water, and exchangeable hydrogen. 
For example, studies of quail fed deuterium-enriched showed that 10 % of all hydro¬ 
gen was exchangeable, 30 % came from ingested water, and the remaining 60 % 
came from food. The overall finding for animals is that 2 H values are primarily 
controlled by diet, which in turn is strongly correlated with 2 H values of local water. 

The reasons why hydrogen isotope tracers work well for studies of long distance 
migration depends on the fundamental chemistry of isotope fractionation. The water 
cycle of evaporation from oceans and precipitation inland involves isotope fraction¬ 
ations that leave behind the light isotopes. During chemical equilibrium fraction¬ 
ation, the vapour phase is enriched in the light isotopes, leaving the liquid phase 
heavier by difference or mass balance. The isotopically heavier liquid phase depos¬ 
its as rain or snow, so that residual cloud-borne water moving inland or upwards is 
isotopically lighter and has lower 2 H values. Larger fractionations accompanying 
cold conditions in Polar Regions magnify some of these patterns, creating low 2 H 
values nearer the poles. This leads to continental-level patterns in the isotope com¬ 
positions of water, in effect a giant isotope map created by the water cycle. Bird 
migrations take place across this chemical landscape, with birds at high latitudes 
having low 2 H values and birds near the equator having high 2 H values. The isotopes 
in the water provide a basic signal that first labels plants during photosynthesis and 
carbohydrate metabolism, then leads to general labelling of the local food web. The 
end result is that organic matter in materials such as bird feathers will have 2 H val¬ 
ues that reflect the local water 2 H values. Migrating birds typically moult and form 
new feathers at the end of the summer, and feathers retain that late-summer isotope 
chemistry until the next year’s moult. 


2.3 The Stable Isotopes of Water on a Spatial Scale 

As surface ocean waters evaporate into the atmosphere, the clouds formed are iso¬ 
topically depleted in 2 H and 18 0 relative to the ocean, resulting in an air mass that is 
similarly depleted relative to the ocean. In turn, as moisture is condensed from 
clouds during precipitation events, that water is isotopically enriched in 2 H and 18 0 
relative to the cloud, leaving the residual cloud mass isotopically depleted in 2 H and 
18 0 relative to the original cloud mass. The process of differential isotope depletion 
during precipitation results in a predictable pattern of depleted isotope ratios of 
precipitation as cloud masses move inland. Since the hydrogen and oxygen in pre¬ 
cipitation become the primary source of H and O atoms incorporated into carbohy¬ 
drates, proteins and lipids during microbe, plant and animal growth, these isotopes 
have the potential to carry a geographically based piece of information that are 
useful in migration studies. 

The International Atomic Energy Agency (IAEA) and the World Meteorological 
Organization (WMO) has been conducting a worldwide survey of oxygen and 
hydrogen isotope content in precipitation. The objective is to collect basic data on 
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isotope content of precipitation on a global scale to determine temporal and spatial 
variations of environmental isotopes in precipitation to provide isotope data for the 
use of environmental isotopes in hydrological investigations. To this primary objec¬ 
tive two other objectives have been added, providing input data to verify and further 
improve atmospheric circulation models, and to study of climate change. Since 
1961, more than 800 meteorological stations in 101 countries have been collecting 
monthly precipitation samples for the Global Network of Isotopes in Precipitation - 
GNIP, IAEA 2006. Precipitation samples are collected in cooperation with WMO, 
national meteorological services and national authorities. The samples are analyzed 
in the IAEA’s Isotope Hydrology Laboratory in Vienna, and in cooperating labora¬ 
tories. Each contributing laboratory is responsible for the accuracy and precision of 
its own analyses. The GNIP homepage can be accessed at: - http://www-naweb. 
iaea.org/napc/ih/index.html. GNIP data may be used freely, provided the source is 
cited as follows: - IAEA/WMO (2006), Global Network of Isotopes in Precipitation, 
The GNIP Database, (http://www.iaea.org/water). From these and other observa¬ 
tions, both temporal and spatial patterns emerged that established the basis for a 
global representation of the distribution of isotopes in water on a global basis. 
Analyses of the spatial distributions of hydrogen isotopes of waters across North 
America and Europe reveal substantial variations in stable isotope ratios, making it 
possible to distinguish precipitation in many geographical locations on the basis of 
their 2 H and 18 0 values. There are no unique stable isotope ratio values for waters in 
a specific geographic location, but rather gradients or bands of different isotope 
ratio values allowing one to distinguish between locations if they were sufficiently 
far apart from each other. Bowen and colleagues extrapolated from this location- 
specific data of stable isotopes in water to construct spatial maps of the predicted 
isotopic composition of water throughout the world so that it is possible to estimate 
reliably the globally averaged 2 H and 18 0 values of precipitation for different lati¬ 
tudes and longitudes using a calculator available at http://waterisotopes.org. 

There are a number of programs available online that can be used to estimate the 
mean annual and monthly H and O isotope composition of precipitation in different 
locations. They are intended to facilitate the use of stable isotope data, amongst 
other uses, in ecological studies, and to standardize the interpolation of precipitation 
stable isotope data. The program for calculating isotopes - The Online Isotopes in 
Precipitation Calculator, OPIC - is based on the work of Bowen and Wilkinson 
(2002), Bowen and Revenaugh (2003) and Bowen et al. (2005). These papers pro¬ 
vide a description of the methods. The data that were used to develop the database 
used by the OIPC are derived primarily from the International Atomic Energy 
Association/(WMO) Global Network for Isotopes in Precipitation (GNIP). To use 
the database it is only necessary to input basic information on location and elevation 
to acquire 2 H and 18 0 values for the specified location. The program is periodically 
updated as new data, including that provided from sources other than the Global 
Network for Isotopes in Precipitation, are added to the database. The OPIC can be 
accessed from the website managed by Purdue University, West Lafayette, Indiana, 
USA, (http://wateriso.eas.purdue.edu/waterisotopes/pages/data_access/oipc.html). 
When the database is used in publications, reference should be made to the data- 
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base, The Online Isotopes in Precipitation Calculator, version XX (http://www. 
waterisotopes.org) and to publications by Bowen and Revenaugh (2003) and Bowen 
et al. (2005). 

Another program that allows the user to download their own isoscape data is 
isoMAP, also based at Purdue University, that allows spatial analysis, modeling and 
prediction of stable isotope ratio variation in the natural environment It comprises 
web-based GIS and software tools that enable users to develop, and implement 
models for isotope distributions. This program can be accessed from http://isomap. 
rcac .purdue. edu: 8080/gridsphere/gridsphere . 


2.4 Deriving Isoscapes in the Absence of GNIP Data 

Numerous studies have shown that precipitation isoscapes drive 5D and 5 18 0 pat¬ 
terns in surface waters and in terrestrial food webs. While the GNIP dataset has 
provided the fundamental spatiotemporal foundation at a global geospatial scale, 
the GNIP stations are often spatially deficient for many regions of the planet that are 
of interest to ecologists. In North America, for example, Mexico has only two GNIP 
stations, which results in a spatial deficiency of water isotope data for that country. 
This spatial deficiency is of concern to researchers interested in large scale Mexican 
isotope hydrology, and to scientists interested in tracking migratory species to and 
from Mexico by using stable isotopes. Overcoming this GNIP spatial limitation 
would require adding spatially dense, long-term precipitation isotope collection sta¬ 
tions across Mexico in order to obtain sufficient isotopic coverage of this topo¬ 
graphically diverse country. This option is not viable for logistical and cost reasons, 
and so other approaches are needed to fill crucial information gaps. A possible 
approach is the application of the generalized globally predictive models for pre¬ 
cipitation isoscapes (www.waterisotopes.org). However, these global models also 
lack sufficient or detailed on-the ground validation to determine how close they 
match or digress from reality at regional or country-wide scales (Bowen and West 
2008). 

In order to overcome this problem Hobson et al. (2009) developed a predictive 
general linear model (GLM) for hydrogen and oxygen isotopic spatial patterns in 
Mexican groundwater and then compared the results to a validation subset of field 
data, as well external data reported in the literature. The GLM used elevation, 
latitude, drainage basin (Atlantic vs. Pacific), and rainfall as the most relevant pre¬ 
dictive variables. The GLM explained 81 % of the overall isotopic variance observed 
in groundwater, 68 % of the variance within this validation subset, and 77 % of the 
variance in the external data set. This predictive GLM is sufficiently accurate to 
allow for future ecological, hydrological and forensic isoscape applications in 
Mexico, and may be an approach that is applicable to other countries and regions 
where GNIP stations are lacking. They hypothesized that the stable isotopic compo¬ 
sition of shallow phreatic groundwater in Mexico might serve as a useful proxy for 
integrating longer term (e.g. 5-10 year) precipitation infiltration inputs (Clark and 
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Fritz 1997). Groundwater infiltration is well-known to be a multi-year integrator of 
seasonally weighted rainfall events (Clark and Fritz 1997; Darling and Bath 1988; 
Rozanski 1985), and hence its isotopic composition closely reflects that of season¬ 
ally weighted long-term average inputs even in some of the most arid regions of the 
world (IAEA 2007). Groundwater was sampled from 234 sites at ~50 km latitudinal 
spacing to obtain high spatial resolution and country-wide coverage for the con¬ 
struction of a groundwater isoscape. Shallow groundwater infiltration in Mexico 
appeared largely unaffected by evaporation and reflected seasonally weighted pre¬ 
cipitation inputs. 


2.5 Use of Stable Isotopes for Migration Studies 

The stable isotopes of H, O, N, C and S are routinely studied in bird migration 
because they are recorders of the dietary sources of the birds and these sources can 
be interpolated or specifically linked to ground proven or large-scale patterns in the 
global landscape and water environments (Hobson and Wassenaar 1997). Depending 
on whether the tissue of the migrating animal is biochemically fixed (feather, hair) 
or dynamic (blood, muscle) these stable isotopes provide fundamental information 
about where the animal has been and what it has been eating. The technology is 
relatively recent in inception and although it has been widely used approaches to its 
use are evolving, techniques are improving and methods of interpretation are being 
updated constantly. Stable isotopes as intrinsic markers do not require a recapture 
of the same animal to obtain data since the required spatial information is retained 
in the animals’ tissues as a result of its stay in a particular habitat. Stable isotopes 
depend on several key requisites in order to function as intrinsic markers. Firstly, 
the animal must acquire a stable isotope into its tissues; since the light stable iso¬ 
topes are key building blocks in animal tissues this is readily achieved. Secondly, 
the migrating animal must move between different environments and retain in its 
tissues measureable isotopic differences that can be linked to diet at previous or 
current locations. This is easily met by species that migrate seasonally across dis¬ 
tinct isotopic landscapes such as between higher northern latitudes and lower lati¬ 
tudes in the southern hemisphere. It also requires that isotopic discrimination 
between the different seasonal landscapes and the tissues being measured is consis¬ 
tent or well-known, or has been empirically measured. Such conditions may not 
always be met with and provide considerable challenge to the interpretation of 
results. Any migration study related to HPAI will need to be conducted collabora- 
tively between virologists, biologists and ecologists with experience in bird migra¬ 
tion and scientists with access to a stable isotope laboratory facility. Stable isotope 
analysis is a costly procedure and it is essential that appropriate links are made that 
will ensure the right approaches to determine migratory connectivity are made and 
that a well-designed project that will answer the questions on virus dissemination 
and wild bird migration. 
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The IRMS used most frequently for SIA is the continuous flow isotope ratio 
mass spectrometer (CF-IRMS). These machines enable relatively low cost analysis, 
a high degree of automation and a high throughput. A technique known as compara¬ 
tive equilibration using CF-IRMS is used for 5D analysis for migration research. 
This approach allows comparative results between different laboratories, ease of 
use and rapid, automated sample throughput of samples (Wassenaar and Hobson 
2003). It requires the inclusion of pre-calibrated keratin working standards along 
with the unknown tissue samples. The hydrogen isotope exchange between ambient 
moisture in the air and the keratin standard and the samples is identical. The com¬ 
paratively equilibrated samples and standards are then isolated from the atmosphere 
and analyzed together in a single session. Measurement of stable hydrogen in tis¬ 
sues is performed on H 2 derived from high-temperature flash pyrolysis. Pure H 2 is 
used as the sample analysis gas and the isotopic reference gas. A high temperature 
Elemental Analyzer and autosampler pyrolyzes samples to H 2 gas. The pyrolysis 
column consists of a ceramic tube filled with glassy carbon chips held at 1250-1350 
°C, followed by molecular sieve gas chromatography column at 80-100 °C. The 
IRMS makes measurements in the following way: gases enter a white-hot region 
where electrons are boiled off and the sample gases are ionized. These ionized mol¬ 
ecules have a positive charge and they are passed through a magnetic field that sepa¬ 
rates them according to their atomic mass and isotopes with the resulting ion beams 
focused into collectors for counting. Computers then calculate the final isotope val¬ 
ues. The stable isotope ratios are then reported as a series of delta (5) notations of 
positive or negative numbers expressed as parts per thousand (%c) relative differ¬ 
ences to an international standard. All 5D results are expressed in units per mil (%c) 
relative to the Vienna Standard Mean Ocean Water using the previously calibrated 
keratin standard. 

One of the most important practical issues in using 5D measurements to track 
large-scale terrestrial animal movement is the problem of hydrogen exchange in 
organic samples It is well established that 5D measurements of feathers (and other 
organic materials) are problematical compared to 5 13 C and 5 15 N analyses due to the 
problem of uncontrolled hydrogen isotopic exchange between “exchangeable” 
organic hydrogen in the feathers and isotopically variable ambient moisture in the 
laboratory environment. If left uncorrected, 5D measurements of the total hydrogen 
in an identical feather will yield different results at analytical laboratories located at 
different geographic locations, as well as over time due to geographical and sea¬ 
sonal changes in the hydrogen isotopic composition of ambient moisture. This 
would make it impossible to compare 5D results between different laboratories. 
Wassenaar and Hobson (2003) showed how keratin 5D analyses by CF-IRMS tech¬ 
nology and the use of keratin standards and comparative equilibration can be used 
to provide cost effective and rapid 5D analyses for migration research to provide 
cost effective and comparable 5D results on feathers and other keratin-based sam¬ 
ples among different laboratories. Three in-house keratin standards were prepared 
from chicken feather cow hoof and bowhead whale baleen. Each keratin type came 
from a single geographic location. All of them were solvent cleaned, dried, cryo- 
genically ground, and homogenized. They were then analyzed 6 times each using 
the offline equilibration dual-inlet method described by Wassenaar and Hobson (2000). 
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The target weight of feather (and in-house keratin standard) required for routine 
online hydrogen isotope analysis by CF-IRMS in our laboratory was 350 ± 10 pg. 
Feather samples were cut from the same position on each feather and weighed using 
a microbalance, and transferred into 4.0 mm 3.2 mm isotope grade silver capsules. 
Capsules containing the keratin standard powder or feather were folded into tiny 
balls. All weights were recorded and the samples were stored in 96-well ELISA 
plates, loosely covered with the lid. Samples and keratin standards were then 
allowed to “air equilibrate” on the shelf with ambient lab air moisture at room tem¬ 
perature for >96 h prior to stable hydrogen isotope analysis. The hydrogen exchange 
process could be sped up by holding the samples and working standards in an oven 
at higher temperatures since showed that complete exchange with moisture occurs 
in less than 30 min at temperatures above 100 °C. For practical purposes, however, 
a 96-h equilibration at room temperature was deemed sufficient and easier, and 
could be easily implemented as a standard laboratory analysis procedure. Following 
96-h comparative equilibration, all samples and standards were immediately loaded 
into the auto sampler carousel of the CF-IRMS and analyzed for 5D. 

The uncorrected hydrogen isotope values of the keratin standards from each 
CF-IRMS auto run were subjected to a least squares regression to derive a correc¬ 
tion formula to be applied to all the keratin standards and feather samples within 
that auto run. A typical run included an in-house keratin standard at the beginning, 
2 keratin standards after every 10-12 unknowns, and 2 or 3 keratin standards at the 
end of the auto run. The authors recommend that researchers performing 6D analy¬ 
ses adopt this protocol in order to provide findings that are comparable among labo¬ 
ratories. They also suggest that the development and distribution of feather and 
keratin working standards that have a wide range of 5D can be done on an individ¬ 
ual basis, or preferably, developed and distributed as a future collaborative effort 
among avian researchers and stable isotope laboratories that provide feather 6D 
analyses. Since CF-IRMS enables 5D analyses to be conducted on small amounts 
of tissue samples it is possible intra-sample variance in the hydrogen isotopic due to 
any biological heterogeneities could exceed interpretations of environmental iso¬ 
tope. To help resolve this, feathers were obtained from captive birds to examine 
isotopic variance expected due to sample size, location, and heterogeneity factors, 
and from selected wild birds to examine isotopic variance due to these and to addi¬ 
tional dietary or location changes during feather growth (Wassenaar and Hobson 
2006). Captive bird feathers were sub-sampled along the vane on either side of a 
single feather at masses of 0.25, 0.35, 0.45, 0.6, 1.0 and 2.0 mg, and along the 
rachis. The results showed consistency of feather 6D measurements across a wide 
range of sample masses. Within-feather 5D isotopic variance for captive and some 
wild birds was as low as ±3 %c for vane material, which corresponds to a geospatial 
resolution of about 1° of latitude in central North America. Intra-sample variance 
for the rachis was ±5%c, with lower 5D values for both wild and captive birds. 
Nevertheless, the authors recommend researchers assess the degree of intra- and 
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inter-sample hydrogen isotopic variation in the selected tissue growth period for the 
species of interest before geospatial interpretations of origin are attempted. 


2.6 Approaches for Determining Migratory Connectivity 

Patterns in stable isotopes derived from the environment are translated into food 
webs and the discrimination between stable isotopes and dietary sources and the 
feather is predictable and constant in time and space. The isotopic values of an ani¬ 
mal are assumed to be representative of the location at which the feathers was grown 
and when the bird moves elsewhere it can be sampled to infer its previous geo¬ 
graphic location. In order to infer the location most likely associated with the feather 
isotopic values it is necessary to calibrate the geographic model using tissues from 
birds of known origin. The accuracy of the IRMS measurements is best when using 
standards with chemical composition similar to the unknown samples (e.g. keratin 
standards to calibrate feather samples). Also, for geographical assignment it is best 
that the standards have similar attributes, i.e. the same species and isotopic values 
as the unknown samples. A calibration dataset should include samples from all 
areas from which the migrant of interest could have originated. However, sampling 
tissues from all potential places of origin is costly and logistically difficult if not 
impossible if the areas are remote and beyond the reach of sampling and collection 
networks. 

The simplest method of assigning a bird to an isoscape is to define the geo¬ 
graphic gradients based on maps of isotopic values, measure an individual and then 
assign it to the area that corresponds to the isotope values obtained. The spatial pat¬ 
terns are generated from indirect sources of information such as rainfall. The earli¬ 
est studies on migration utilized this technique (Hobson and Wassenaar 1997) to 
use 5D to study migratory birds and it is still in common use (Hobson et al. 2006, 
2007). The method is simple to understand and apply and no additional computation 
is required to place birds in their geographical locations. More advanced methods 
for determining the origin of migratory animals require some degree of statistical 
manipulation. Assignment methods require the definition of all possible locations of 
origin, then characterizing these locations with the distribution of stable isotope 
data. The characterizing should, if possible, be derived from measurements of birds 
known to have grown tissues at those locations and all possible regions of origin are 
sampled. Once a calibrated assignment model is obtained, stable isotope measure¬ 
ments from individuals of unknown origin are used to determine their origins. 

Although migratory individuals have only to be captured once in order to deter¬ 
mine their migratory movements the assignment of their origin is challenging and it 
is necessary to make inferences about the history of an individual using the stable 
isotope measurement. Table 2.5 details some of the different modelling approaches 
that have been applied to estimate origin when using stable isotopes. 


Table 2.5 Methods for assigning animals of unknown origin using stable isotopes 
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References 

Hobson et al. 

(2007) 

Kelly et al. (2002) 

Hebert and 

Wassenaar (2005) 

Kelly et al. (2005) 

Incorporates 
sources of error 

No 

No 

No 

Some 

Disadvantages 

Does not 

incorporate spatial 
variability of 
isotopes 

Does not 

incorporate spatial 
variability of 
isotopes within 
regions 

Does not 
incorporate most 
error, no a priori 
hierarchy for 
multiple isotopes 
does not provide 
degree of certainty 
for branching 

decisions 

Does not 

incorporate most 

error, simply 

assigns regions for 

individual based on 

highest likelihood 

value; requires 

sampling all 

potential regions of 

origin 

Advantages 

Easy to implement 

Easy to implement 

Can be applied to 
multiple isotopes. 

Does not require 

distributional 

assumptions 

Can be applied to 
multiple isotopes. 
Provides probability 
of assignment for a 
given individual, 
easy to implement 

Description 

Isotope value of 
animal assigned to 
area based on 
isoscapes on map 

Origin inferred on 
regression of isotopes 
on latitude and/or 
longitude 

Origin inferred on 
basis of hierarchical, 
discrimination based 
decision rules 

Origin inferred from 
probability density 
functions for isotope 
values from given 
regions 

Model type 

Map lookup 

Linear regression 

Classification trees 

Likelihood-based 

assignment 
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References 

Norris et al. (2006) 

Wunder and Norris 

(2008) 

Wunder (2010) 

Incorporates 
sources of error 

Some 

Yes 

Yes 

Disadvantages 

Does not 
incorporate most 
error. Simply 
assigns region for 
individual based on 
highest posterior 
probability; requires 
sampling all 
potential regions of 
origin 

Requires intensive 
computing 

Requires intensive 
computing 

Advantages 

Can be applied to 
multiple isotopes 
and can utilize 
non-isotopic 
information 

Can be applied to 
multiple isotopes, 
incorporates multiple 
sources of error and 
provides a range of 
assignments for a 
given individual 

As above, 

incorporates known 
variance sources 

Description 

Same as above but 
adds prior information 
on isoscapes using 
Bayes’ Rule 

Same as above but 
adds sources of error 
associated with 
isotope data 

Models stochastic 
error process over 
mean field surface 

Model type 

Likelihood with 
priors 

Stochastic extension 
of likelihood 

Probability Surfaces 
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2.6.1 Classification Trees to Predict Origins 

Ecological data are often complex, unbalanced, and contain missing values. 
Relationships between variables may be strongly nonlinear and involve high-order 
interactions. The commonly used exploratory and statistical modelling techniques 
often fail to find meaningful ecological patterns from such data. Classification trees 
are statistical techniques ideally suited for exploring and modelling such data. They 
have a number of advantages over discriminant function analysis and linear regres¬ 
sion, both of which are often used in stable isotope assignment studies. Classification 
trees represent a statistical technique well suited for modelling ecological data since 
the model development is hierarchical and based on logical if-then conditions and 
are both nonparametric and nonlinear (De’ath and Fabricius 2000). Hebert and 
Wassenaar (2005) used the classification tree approach to determine if ducks origi¬ 
nating from different geographic areas had unambiguous multi-isotopic signatures. 
They used a multi-stable isotope analysis (5 34 S, 5 2 H, 5 13 C, 5 15 N) of secondary feath¬ 
ers from wild duck (Anas platyrhynchos ) and northern pintail (A. acuta) ducklings 
from 52 sites in western North America. Ducklings from Alaska, northern Canada, 
the Prairies and California could be distinguished based upon their feather isotope 
values. Classification trees were developed by repeatedly splitting the data via algo¬ 
rithms that partition the data into mutually exclusive groups (De’ath and Fabricius 
2000). Geographic patterns in feather isotopes were related to natural gradients pro¬ 
duced by biogeochemical cycles and anthropogenic factors such as agrochemical 
usage. 


2.6.2 Likelihood-Based Methods to Predict Origins 

Discriminant analysis can be used to classify a sample into two or more classes 
(regions) and has been employed in determining the origin of birds. Parameters for 
region-specific likelihood functions are estimated from isotope data collected in 
each predefined region. The likelihood functions for each region are then evaluated 
for the isotope value measured for an individual of unknown origin and that indi¬ 
vidual is placed in the region with the highest-valued likelihood (or probability). 

In studies on golden plovers (Prosser et al. 2006) interpolated maps were used of 
annual average 5D values collected at stations across the wintering grounds of each 
of two species of plover to generate a priori predictions for changes in 5D in feathers 
between breeding and wintering grounds (IAEA 2001; Bowen et al. 2005). Because 
there was no information on the specific wintering location of each bird, the expected 
range of 5D feather values was evaluated from these maps (Bowen et al. 2005). The 
annual averages of deuterium in the precipitation across wintering areas were used 
to predict feather 5D values. To assign individuals to group of origin, the principle 
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of “assignment with exclusion,” was used, a method adopted from population genet¬ 
ics that uses an exclusion criterion to reject unlikely sources. The cross-validation 
procedure in discriminant analyses as an assignment method was used to classify 
feathers into known group of origin based on the three isotopic values. The limited 
number of sample contributed to low power of some assignments, but power was 
increased by using three isotope values for each sample. 


2.6.3 Migration Studies Using Stable Isotopes 

The first premise for using stable isotopes is that there are differences in the spatial 
distribution of the isotopes from where the migrants obtained their diet. The stable 
isotopes 5D and 5 18 0 have the most predictable geographic gradients, but gradients 
for other light isotopes are less well known. Geographic patterns of stable isotopes 
are inferred from a number of sampling points across the landscape. The 5D values 
in precipitation have been modelled for North America (Hobson and Wassenaar 
1997) where the relationship between 5D in feathers and growing season precipita¬ 
tion was examined for neotropical migrant songbirds breeding over a continent¬ 
wide isotopic gradient. The 5D values were determined on feathers of 140 
individuals of 6 species of wild insectivorous forest songbirds (Setophaga ruticilla, 
Empidonax minimus, Vermivora peregrinus, Catharus ustulatus, Seiurus aurocapil- 
lus, Hylocichla mustelina) taken from 14 breeding locations across North America. 
The 5D of feathers was strongly correlated with the 5D of growing season precipi¬ 
tation at breeding sites across North America. As feather hydrogen is metabolically 
inert after growth, this relationship was then used to assess the breeding origins of 
wintering migrants. Deuterium values of feathers from 64 individuals representing 
5 species of migrants ( Helmitheros vermivorus, Wilsonia citrina, Hylocichla muste¬ 
lina, Dumetella carolinensis, Seirus aurocapillus ) at a wintering site in Guatemala 
were consistent with those predicted from the known breeding ranges of these spe¬ 
cies. Bowen et al. (2005) described a method for interpolation of precipitation iso¬ 
tope values to create global base maps of growing-season (GS) and mean annual 
(MA) 5D and 5 18 0. The use of these maps for forensic application was demon¬ 
strated using previously published isotope data for bird feathers in North America 
and Europe. The precipitation maps show that the greatest potential for applying 
both hydrogen and oxygen isotope forensics exists in mid- to high-latitude conti¬ 
nental regions, where strong spatial isotope gradients exist. They showed that 
feather 5D/precipitation 5D relationships have significant predictive power both in 
North America and Europe, and zones of confidence for the assignment of origin 
could be described using these predictive relationships. These maps are available in 
GIS format at http://www.waterisotopes.org i 
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2.6.4 Determining Migratory Connectivity for Waterfowl 
in Asia 

Although Asia is a prime focus for HPAI, and despite considerable interest in con¬ 
servation and disease surveillance issues involving waterbirds in that region inves¬ 
tigations on migratory connectivity in wild birds that might be implicated in the 
spread of the disease are few (Chang et al. 2008; Perez et al. 2010). The great cor¬ 
morant ( Phalacrocorax carbo sinensis) is a potential carrier of HPAI to Taiwan and 
to better understand the migration of this bird, SIA studies were conducted on sev¬ 
eral populations from different breeding sites and overwintering grounds (Chang 
et al. 2008). Rather than capturing and removing feathers from live adult birds, 
newly dropped feathers found lying on the ground including primary, secondary 
and tail feathers at collected from three breeding sites in China (Qinghai Lake), and 
Russia (Ussuri River, and Chita Peninsula) and overwintering grounds in Taiwan 
(Kinmen Lake). All the collected feathers were newly dropped because these feath¬ 
ers were neither weathered, nor damaged, nor covered by a layer of dirt or dust. 
They were collected from different individuals by collecting them under widely 
separated nests at breeding sites and under separated roosting trees. Although 5D 
values of feathers have typically been compared with the 5D of precipitation to 
infer migration, the distributions of 5D and 5 18 0 in precipitation across Asia are not 
as precisely known as in North America and Europe as there are fewer data collec¬ 
tion sites and the frequency of data collection is more sporadic in Asia. This makes 
it more difficult to infer the origins of birds by comparing the 5D values of precipi¬ 
tation and of feathers in Asia alone hence, a direct comparison of 5 13 C, 5 15 N, 5 18 0 
and 5D values of feathers collected at both wintering and breeding sites was used to 
infer breeding locations of P. c. sinensis wintering at Kinmen. Analysis showed that 
it was unlikely that P.c. sinensis wintering at Kinmen came from the breeding sites 
in China and Russia. 

Analysis of the data collected from the P.c.sinensis wintering at Kinmen, pro¬ 
duced nine clusters plus two outliers and approximately 65 % of the birds overwin¬ 
tering at Kinmen come from four breeding populations. Using 5 13 C, 5D, 5 15 N or 
5 18 0 values, excluded the possibility that cormorants wintering at Kinmen came 
from breeding populations at Qinghai Lake, Ussuri River or the Chita Peninsula. 
The general locations of possible breeding sites were based on the comparisons of 
feather 5D values with the mean 5D of precipitation during the breeding season 
according to the best available 5D isocline maps of Asia, from which the data are 
spatially and temporally discontinuous. This comparison suggested that approxi¬ 
mately 50 % of the P.c.sinensis wintering at Kinmen come from five breeding popu¬ 
lations from the region around Lake Baikal or an area encompassing the Amur, 
Khabarovsk and Primorsky regions of Russia. The study showed that SIA has the 
potential to provide migratory information at the population level, to uncover the 
possible existence of previously unknown breeding populations, and to make pre¬ 
dict migratory and possible future transmission pathways of HPAI. Moreover, com¬ 
bining multivariate analyses of stable isotopes with analyses of the phylogenic trees 
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of birds and HPAI may provide an opportunity to better understanding the interac¬ 
tion of bird migrations and the spread of HPAI. 

In another study (Perez et al. 2010) feather sample were collected from Bar¬ 
headed Geese ( Anser indicus ), Whooper Swans ( Cygnus cygnus ), Mongolian Gulls 
(Larus vegae mongolicus ), Curlew Sandpipers ( Calidris ferruginea) and Pacific 
Golden Plover (Pluvialis fulva) at seven sites in northern Mongolia. After cleaning 
to remove oils, stable-hydrogen isotope (5D) analysis of one feather per individual 
was conducted at the stable isotope facility at the National Hydrology Research 
Centre, Saskatoon, Saskatchewan, Canada. Isotope analyses followed the compara¬ 
tive equilibration technique described by Wassenaar and Hobson (2003). Deuterium 
isotope ratios were expressed in 5 notation in parts per thousand (%c) relative to the 
Vienna Standard Mean Ocean Water Standard Light Antarctic Precipitation 
(VSMOW-SLAP) scale. Measurement precision based on replicate measurements 
of within-run keratin standards was estimated to be of the order of ±2 %c. Estimates 
of continental deuterium patterns in mean growing-season precipitation (5Dp) for 
Asia were derived from Bowen 2009 (www.waterisotopes.org). Continental 5Dp 
data show an expected general decrease in deuterium abundance in precipitation 
with increasing latitude and altitude in Asia and form the basis for predicting the 
expected abundance of deuterium in feathers (5Df) from the average 5Dp (Hobson 
and Wassenaar 1997; Hobson 2008). Three major feather isotopic source areas in 
Asia were identified: Region A - northern Asia (<-138%c), B - catchment area 
(-138 to -104%c) and C - southern Asia (>-104%c) were delineated. The isotopic 
interval of 34%c presented for Region B corresponded to latitudes ranging from 
approximately northern China to mid-Russia, and Regions A and C were north and 
south of those latitudes, respectively. These arbitrary cut-off values in 5Df were 
chosen because they bracketed Region B between the lowest and highest 95 % Cl 
out of the expected mean annual 5Dp values for the latitudes, longitudes and alti¬ 
tudes of the sampling sites. Once more precise information on the relationship 
between 5Dp and 5Df can be established for Asia, more refined assignment will be 
possible (Hobson et al. 2009). Feathers from adult Pacific Golden Plover 5Df had 
values of-33%o, clearly indicated that these feathers were moulted in southernmost 
Asian latitudes, suggesting that the feathers were probably grown on or near the 
wintering grounds prior to spring migration. Locally grown feathers from adult 
Barheaded Geese and Whooper Swans had more depleted 5Df values than any other 
group because they live in a complex of lakes and ponds fed by high elevation run¬ 
off, and the waters are more depleted in deuterium than adjacent areas. Curlew 
Sandpiper had 5Df values consistent with Region B indicating isotopic equilibrium 
with diet on their sampling area. The study demonstrated that stable isotopes could 
be used to improve knowledge of migratory and moulting patterns of wild water- 
birds in Asia in a rapid and cost effective manner. However, additional studies are 
required since the 5Dp isoscape for Asia is so poorly described. Ground-truthing 
the relationship between 5Df and 5Dp over a large geographic gradient in Asia will 
provide more information to derive appropriate isotopic calibrations (Wunder and 
Norris 2008; Wunder 2010). 
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Arrival time on breeding or non-breeding areas is likely to be of interest in epi¬ 
demiological studies exploring dissemination of HPAI by migratory birds. Precisely 
assessing the arrival time of individuals can be difficult, but by measuring carbon 
stable isotope turnover in avian blood it is possible to estimate arrival time for birds 
switching from one habitat to another (Oppel and Powell 2010). Stable carbon iso¬ 
tope ratios (5 13 C) in blood assimilate to a new equilibrium following a diet switch 
according to an exponential decay function. Stable carbon isotope ratios in tissues 
reflect the isotope ratios of food sources (Hobson and Clark 1992). When migratory 
birds switch habitat their new diet will have a different isotopic signature if the two 
habitats are isotopically distinct (Peterson and Fry 1987). This switch to a new diet 
causes the isotope ratios in blood to change gradually over time until they reflect the 
isotope ratio of the new diet (Hobson and Clark 1992; Evans Ogden et al. 2004; 
Morrison and Hobson 2004). The rate of change in 5 13 C is tissue dependent. For 
example, blood plasma generally assimilates to the stable isotope ratio of a new diet 
within a few days (Hobson and Clark 1992), whereas whole blood or the cellular 
fraction of blood (red blood cells, RBC) turn over within several weeks (Bearhop 
et al. 2002; Evans Ogden et al. 2004; Morrison and Hobson 2004). Several experi¬ 
mental studies have determined that isotopic turnover in blood closely follows an 
exponential decay function (Evans Ogden et al. 2004; Carleton and Martinez del 
Rio 2005). This relationship can be used to determine the time an animal switched 
diets if the isotope ratios of the old and the new diet are known. 

A study of the arrival time of Eider Ducks (Somateria spectabilis) at breeding 
grounds (Oppel and Powell 2010) was estimated using only a single tissue (Phillips 
and Eldridge 2006). Since the turnover rate of stable carbon is not known in many 
bird species a mass dependent turnover rate constant was utilized in calculating the 
arrival time (Carleton and Martinez del Rio 2005). Data from whole blood and RBC 
measured in captive experiments were used to validate the approach of estimating 
the time since a diet switch. In each experiment, birds were switched from one iso¬ 
topically distinct diet to another differing by at least 3 % in 5 13 C (Hobson and Clark 
1992, 1993; Bearhop et al. 2002; Evans Ogden et al. 2004). This provided data to 
develop a formula that could be used to determine arrival time. If the isotope signa¬ 
tures of a tissue from both the old and new environment are known and birds are 
captured at an unknown time after arrival, then blood 5 13 C can be used to determine 
the time since elapsed since arrival. The isotopic signature of each diet can be char¬ 
acterized by sampling blood from birds that have been feeding in either environ¬ 
ment for several weeks, or by using tissues that are metabolically inert after growth. 
Feathers sampled from the bird may be useful to determine 5 13 C of the previous 
environment if the moulting strategy of the bird is sufficiently known and certain 
feathers are always grown in the previous environment. Plasma or RBC sampled 
from a bird captured at a later stage can be used to determine new 5 13 C once birds 
have reached equilibrium with the new diet (Morrison and Hobson 2004). For 
breeding grounds, eggshell membranes collected from old nests provides a simple 
way to characterize the diet of birds (Oppel and Powell 2010). 
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Stable isotopes may also provide clues to help locate nonbreeding populations of 
birds that have unknown winter ranges. For instance, no records exist for the Coastal 
Plain Swamp Sparrow (Melospiza georgiana nigrescens) from October to late April 
because of the difficulty in tracking individuals between seasons. Stable isotope 
analyses of C, N, and H were used to predict where moult occurs and then those 
areas were then searched for individuals of this species (Greenberg et al. 2007). 
Feathers were clipped from the crown and rump and lower back of the birds. 
Analysis of the crown feathers was used to predict the 5D value at the unknown site 
of winter migration. Selection of the sites was based on the estimates of isotopes in 
precipitation reported by Bowen and colleagues (www.waterisotopes.org). The 
5 13 C and 5 15 N of rump feathers were consistent with the Sparrow moulting in more 
saline marshes. The values for the same isotopes from crown feathers revealed that 
winter moult in the Sparrow probably occurred in similar coastal brackish habitats. 
The 5D of crown feathers indicated that pre breeding molt occurred at latitudes 
between South Carolina and Virginia. A subsequent search of this region located 
specimens of the birds, all of which were found in North Carolina or southeastern 
Virginia. The Sparrows were found predominantly in brackish marshes similar to 
their breeding habitat. On the basis of these observations, it appears that Coastal 
Plain Swamp Sparrows undergo a short southerly migration to a coastal region with 
substantially warmer winter conditions. This study was the first to make a specific 
geographic prediction based on stable-isotope analysis and to test the prediction in 
the field. 

Many birds that pass through Africa on their migration cross the Sahara or 
Arabian deserts without delay, then stopover in northern tropical areas for several 
months before resuming their flights to the south. During the stopover they com¬ 
plete a partial or complete moult (Yohannes et al. 2007). Stable nitrogen (5 15 N), 
carbon (5 13 C) and hydrogen (5D) isotope profiles in feathers of nine migratory bird 
species trapped in Kenya were examined to test the extent to which they were seg¬ 
regated, geographically or by habitat, during their earlier autumn migration stop¬ 
over in northeast Africa. The aim was to determine if isotopic differences between 
species varied between years, and whether the isotope profiles of individual species 
appeared to be consistent. The relationship between mean feather 5 13 C, 5 15 N and 
5D assorted the migrants into several clustered groups. Similar feather isotope val¬ 
ues among successive years revealed that each species tended to return to the same 
or similar stopover areas and selected habitat and diet that generated similar isotopic 
signatures. However, the stopover sites that the birds used were not identified. 
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Chapter 3 

Practical Considerations 


3.1 Sample Collection and Tissue Preparation 

In assessing the value of stable isotopes for tracing animal migration the O and H 
isotopes in the environment are “global spatial” in their resolution as the isotope are 
related to hydrological and meteorological processes that are seasonally and spa¬ 
tially predictable over many years on a regional, continental and global scale. This 
enables interpolation into regions where no long term data are in existence (Bowen 
et al. 2005). Isotopes of C, N, and S are local spatial in character as they do not vary 
over the landscape as H and O. investigations into large scale migrations is therefore 
most fruitful using H and O isotopes, while the resolving power at the local habitat 
level can be improved with data on the local spatial isotopes. 

Two types of samples that can be used for analysis; the first type is fixed tissue 
consisting of body material that once formed is metabolically and therefore by 
extension, isotopically inert. This material will give an indication of the isotopic 
composition pertaining in the local environment at the time that the tissue was 
formed. Tissues of this type include keratinous material of claws, nails, hair and 
feathers. Bird feathers are particularly well-suited for use in tracking migration as 
the fully formed feather vanes and rachis do not change chemically or isotopically 
as the bird moves away from the site of feather formation, on its migration pathway. 
The feather is often grown over a short period of time at a specific site such as a 
breeding site or wintering ground and so will retain the isotopes present in the diet 
at that location. It is essential that the researcher has adequate knowledge of the 
ecology and physiology of the birds being studied. For instance, slowly growing 
feathers or claws that grow as an animal migrates, ingesting different diets en-route 
means that a growing feather will record a variation in isotope as it progresses. In 
bald eagles where flight feathers grow as they migrate from the northern to the 
southern USA, there is intrasample variation where the feather show a negative 5H 
value in the oldest part of the feather at the wing tip, while at the base of the feather 
the 5H values are more positive. In this case the feather is not equivalent to a fixed 


© IAEA 2016 

G.J. Viljoen et al., Stable Isotopes to Trace Migratory Birds and to Identify 
Harmful Diseases, DOI 10.1007/978-3-319-28298-5_3 


35 



36 


3 Practical Considerations 


tissue and provides instead a dynamic spatial recorder. In the case of a migratory 
bird in which the feather is formed at the natal location there would be no significant 
intrasample differences in 5H that would indicate where the bird had moved, even 
though captured far from the region in which the feather was formed. Archived or 
museum material is also a potential source of investigation for stable isotope analy¬ 
sis, but care must be taken in using such material if it has been preserved or con¬ 
taminated in some way whilst in storage. 

The second type of sample material is dynamic tissue, metabolically active tis¬ 
sues in which the stable isotope composition reflects the current, active dietary habi¬ 
tat. These tissues are the blood, muscle, liver etc. There have been few studies using 
such tissues to plot migratory connectivity. Most studies have concentrated on mea¬ 
suring 5 13 C and 5 15 N. The turnover rate for isotopes in different tissues varies from 
a few days in liver and blood, to weeks in muscle to a lifetime in bone collagen. 

An important consideration is whether the sample can be obtained non-lethally 
by plucking or cutting feathers, nails or hair so as not to inconvenience the animal. 
Animal care guidelines should always be exercised when sampling live animal in 
the field. Collecting a sample for analysis is a major undertaking and it is essential 
that due consideration has been given for the animal’s biology and ecology to ensure 
that appropriate samples are collected. 


3.1.1 Questions of Sample Heterogeneity 

For birds, flight feathers or contour feathers can be sampled for isotopes, but to if 
replicate samples are obtained what extent would isotopic values vary from the 
isotopic range that defines the bird’s spatial distribution? There will be individual 
and population level intersample isotopic heterogeneity. There could be minor vari¬ 
ation in values between contour feathers grown by an individual at a single location. 
There would as also be intersample isotopic variation for the same feathers obtained 
from a population of birds that grew the feathers at the same location. This sort of 
heterogeneity is usually much less than the large scale geospatial isotopic patterns. 
Intrasample isotopic heterogeneity is the variance seen in a single discrete sample; 
an example of this is as mentioned above for the bald eagle. 


3.1.2 Preparation of Feather Samples 

Samples collected from the field are likely to be dirty and greasy. Dirt can be 
removed by washing in distilled water but surface oils on feathers will need to be 
removed by organic solvents since the composition of H in the oils differs from that 
contained in the feather keratin, thereby giving false data on isotope composition. 
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Table 3.1 Procedure for preparing feather samples for 5D stable isotope analysis 
Procedure 

1. In order to remove surface oils clean feathers in 2:1 v/v chloroform/methanol by soaking for 
24 h and rinsing twice in same mixture. Air-dry in a fume hood for at least 48 h 

2. Cut off a small amount of feather vein, not the central rachis, for analysis. Always cut 
samples from same region on different feather samples to ensure consistency (e.g. sample at tip). 
Cut feathers using stainless steel scissors and forceps 

3. Clean instruments using methanol and tissue wipes and allow to dry. Do not use acetone 

4. Clean and calibrate microbalance 

5. Tare a silver capsule* handling it only with instruments, remove and set on a metal surface. 
Use the smallest size that will safely hold a sample (3.5 x5.0 mm) 

6. Transfer a small amount of feather to capsule using forceps. Reweigh and add or subtract 
material until the target sample weight of 350 ± 10 pg is achieved. Reference materials must be 
weighed to a comparable elemental mass as the samples 

7. Seal the capsule by crimping the end shut with straight edge forceps and then fold down 
tightly as if rolling down a bag. Use forceps to gently compact the capsule into a small tight 
cube or ball. There should be no loose sides, stray edges or feathers poking out. Flattened 
samples can jam the autosampler 

8. Record sample weight and sample name in an Excel spreadsheet conformed in a 96-well 
ELISA template format. Put capsule in appropriate position in ELISA plate. Clean all 
instruments with methanol and tissue wipes after completing each sample. Put lid on ELISA 
plate and secure with rubber bands or tape to ensure samples do not jump out of wells in later 
handling. Record sample name etc. and weight and its position in the ELISA plate. Transfer 
information to an isotope laboratory submission form 

*Silver Capsules must be used for analysis of 5 18 0 and 5D. They can be obtained from: (1) Silver 
capsules for solid samples, 3.5 x 5 mm, 041066 240-054-00; Costech Analytical Technologies Inc.; 
26074 Avenue Hall, Suite 14; Valencia, CA 91355; USA or (2) D2003 - Silver Capsules Pressed 
3.5 x 5 mm; OEM Parts: E12050, 041066, SC1082; Elemental Microanalysis; 1 Hameldown Road; 
Okehampton, EX20 1UB; UK 


Samples collected for stable isotopic analysis must be stored properly both before 
and after preparation to ensure there is no degradation. Many laboratories can pro¬ 
vide a service but if the researcher wishes to reduce costs then cleaning and weigh¬ 
ing can be undertaken. This is not to be undertaken lightly as accurate weighing to 
±10 pg is required using analytical microbalances. A recommended procedure for 
cleaning and preparing feather samples for analysis is detailed in Table 3.1. 

In order to be successful, the comparative equilibration approach for SIA relies 
on the long term availability of keratin and organic tissue working standards because 
it is essential that there is rigid QA and comparability of results between laborato¬ 
ries to ensure accuracy of results. Although there are some primary, certified organic 
standards - e.g. IAEA-CH-7 which has been checked for isotopic homogeneity and 
has a recommended 5 13 C value determined by international calibration - they do not 
contain exchangeable hydrogen and are not suitable reference material for calibrat¬ 
ing 5D in feathers. Three keratin standards have been prepared at the laboratories of 
Environment Canada, Saskatoon, Saskatchewan, Canada from cryogenically ground 
and isotopically homogenized chicken feathers (CFS), cow hoof (CHS) and bow- 
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head whale baleen (BWB-II). They were calibrated according to Wassenaar and 
Hobson (2003). CFS was obtained from a single batch of 2 kg of chicken feathers; 
CHS was prepared from a hoof of a cow obtained from an abattoir and BWB-II was 
obtained from the University of Alaska. All of the keratins were solvent cleaned 
using methanol/chloroform, air-dried and cryogenically homogenized to about 
0.5 kg. They were passed through <100 pm sieves and further homogenized to 
ensure isotopic homogeneity at the 100 pg level. These standards were sufficient to 
last for several years of use in a single laboratory. The results for 5D were CFS, 
-147±5 %, (VMSOW), CHS, 187±2 % (VMSOW) and BWB-II, -108±4 % 
(VMSOW). 

It is essential that isotopic working standards for 5D in keratinous materials are 
produced to provide for the long-term needs of migration and ecological research. 
Such materials would need to be: (i) in sufficient quantities to allow their use by 
many different laboratories working on stable isotopes for up to 10 years; (ii) be 
certified for 5D homogeneity to <100 pg and (iii) be certified for an isotopic range 
of at least 200 %. This need for standardization will become more important as stud¬ 
ies on stable isotope for migration and other reasons increases to become a main¬ 
stream analytical tool. 


3.2 Other Stable Isotopes for Migration Research 

Stable oxygen isotopes 5 18 0 could also be used for migration research, especially as 
keratins have no exchangeable oxygen so there is no need for comparative equili¬ 
bration procedures. The disadvantage is that the 5 18 0 range for tissues is small 
(approx. 15%c) and the analytical error for measurements are high so that precision 
regarding geospatial information is lost. Carbon (5 13 C) and nitrogen (5 15 N) isotopes 
can also be used in migration studies for local-spatial analysis, to refine populations 
or indicate habitat type. 


3.2.1 Sampling Instructions for Water 8 18 0 and 8 2 H - Rivers/ 
Lakes/Groundwater 

Materials 

Twenty five to sixty milliliter of HDPE (plastic) bottles with caps. Only a few ml of 
water is required for the analysis in the laboratory, but very small bottles are awk¬ 
ward for field staff to clearly write labels on so for logistical reasons use bottles of 
a capacity of 2 ml. Do not use glass bottles since they can shatter in transport. 
Evaporation and leakage during storage must be avoided. This can occur by mishan¬ 
dling or from loose caps. 
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Fig. 3.1 Nomenclature of wing feathers of birds (Sample from the CRP D32030) 


Method 

Fill sample bottle with raw river, ground or lake water to overflowing. Cap tightly. 
Before sampling ensure the sample bottle is clearly and labeled with site or station, 
and dated using permanent ink. Ensure label does not smear. Do not freeze. Store at 
room temperature. No preservation is required. Samples should be shipped ground 
or air with adequate protection (box with foam chips, bubble wrap) and sent to the 
isotope laboratory of choice. 


3.2.2 Sampling Instructions for Bird Feathers for Use in SIA 

Materials 

Steel laboratory tweezers. 

Steel Scissors 

Paper Envelope - Sealable paper envelopes are best to use for storage of bird feath¬ 
ers. Stored feathers should not be wet. Sealable plastic bags can be used as well. 

Long wing feathers: Cut long feather from wing (Fig. 3.1) with pair of scissors 
(preferably PI or P2 feathers). 

Down feathers: Pluck down feathers (5-10) with tweezers from wing, or breast. 

If possible check bird body for Blood keels (feathers in growth) (Fig. 3.2): Pluck 
blood keels (down feathers in growth) with tweezers from breast, or cut one large 
blood keel (large feather in growth) with pair of scissors from wing. Feathers 










40 


3 Practical Considerations 



Fig. 3.2 Blood keel feathers in growth 

should be dry. No visible infestation of feather inhabiting parasites. Put all feath¬ 
ers from one individual bird in the same envelope. Close, seal and label the enve¬ 
lope. Note which samples it contains (large feather, down feathers, blood keel(s); 
bird species, [if possible also: sex, growths status (juvenile, adult)], date and 
location of sampling, name of person sampling. The samples can be stored at 
room temperature. No preservation is required. They should be shipped by 
ground or air with adequate protection (bubble wrap envelope) and sent to the 
isotope laboratory of choice. 


3.2.3 Sampling Wild Birds in the Field 

A comprehensive guide and strategy for sampling wild birds is given by Whitworth 
et al. (2007) and that manual should be referred to in dealing with HPAI (http:// 
www.fao.org/docrep/010/al521e/al521e00.htm). The following section highlights 
a few relevant points. 

On arrival at the site, evaluate the extent of the mortality, if appropriate, includ¬ 
ing number of birds, species directly involved, other wildlife or domestic animals 
involved, and geographic range of mortality. This information should be recorded in 
a Sick or Dead Bird Sample Collection Log. In addition to preparing for animal 
sample collection, also consider collection of other environmental samples includ¬ 
ing water, soil, vegetation or other elements from which it might be possible to 
identify HPAI virus. If it is possible to get a GPS locations that characterize the 
extent of the die-off area, this is preferable to a general verbal description. 

Wear the appropriate level of personal protective equipment, based on the situa¬ 
tion being investigated. Try to minimize direct contact with dead birds and always 
keep animals away from the face, wear vinyl or latex gloves if handling a dead bird. 
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The best method for collecting a dead bird is to invert a plastic bag around your 
gloved hand and then surround the animal with the bag so that the bird is not directly 
touched. Seal the bag tightly (double bag if required for strength and cleanliness) 
and clearly and indelibly label the bag with an Animal Identification Number (which 
must match the number entered in the Sick or Dead Bird Sample Collection Log), 
species, date, time and location. If more than one species has been affected, collect 
several specimens of each for diagnosis. In general, carcasses of birds that have 
been dead for less than 24 h (fresh carcasses) are sufficiently adequate (moribund or 
viraemic birds are best) for diagnostic purposes. In colder climates, carcasses may 
last in relatively good condition for longer periods of time; in warm climates, car¬ 
casses will decompose faster. 

When possible, fresh carcasses should be refrigerated (NOT frozen); a decom¬ 
posing carcass is desiccated, bloated, green, foul smelling and has feathers that pull 
out easily. To increase diagnostic value, fresh carcasses must be transferred to the 
appropriate veterinary or pathology facility and examined as soon as possible. In 
field settings and/or far from appropriate diagnostic facilities, collect samples on 
site and place them in an ice chest or cooler. Keep carcasses away from refrigerators 
used for animal or human food. 


3.2.4 Sampling Strategy for H5N1 

For each affected species, select up to 3 birds that have most recently died (less than 
24 h) or more if practical, up to 3 sick birds (suffering respiratory, neurologic or 
gastrointestinal disease or moribund) and up to 3 apparently healthy birds in direct 
contact with currently sick birds. If possible, also conduct a survey of other live 
birds that share the same habitat (cloacal swabs and/or tracheal swabs only). Priority 
should be given to birds that share wetlands with affected birds since the main mode 
of transmission of AI virus is probably faecal contamination of water, shores or 
banks. It is best to collect as many carcasses as possible and to place them in a cen¬ 
tral location for processing. Removal of dead birds from the site may also help 
prevent secondary contamination of scavengers or the environment. It is very impor¬ 
tant to complete the Sick or Dead Bird Sample Collection Log as carcasses are 
being collected and processed. If possible, try to collect and examine sick animals 
as well as newly dead birds, making sure that you have the appropriate permits to 
capture live samples. If there are too many dead to be able to individually bag and 
label, try to bag or examine well preserved animals that will be most useful for 
diagnostic purposes, and keep these separate from decomposing carcasses. If pos¬ 
sible, transport carcasses (in sealed bags) in a space well separated from the occu¬ 
pants of the vehicle. If you are working in a remote area, you may have to perform 
field necropsies on site. 

Recommended ornithological information to be collected during surveillance 
programmes or field assessment of wild bird mortality events. 
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All birds from which samples are taken should be identified to species. Where 
clearly distinguishable subspecies or discrete populations exist, as for some geese, 
this information should also be collected and reported. Age and sex should be 
recorded wherever possible. Close collaboration with ornithologists in the capture 
and sampling of live birds not only facilitates identification of birds but also gives 
the opportunity to collect additional information on the sampled live birds (such as 
weight, age, sex and condition), which are important to developing a better under¬ 
standing of viral ecology and epidemiology. Standard protocols exist for the collec¬ 
tion of such data through national ringing schemes [details of which are available 
for Europe, for example, via EURING (www.EURING.org)]. Recording individual 
ring numbers in the reporting spreadsheet provides a means of accessing these data 
for future analysis. To provide an audit of identification, it is highly desirable that a 
clear digital photograph is taken of each sampled bird (especially those found dead 
and/or not identified by ornithologists) and stored at least until confirmation of labo¬ 
ratory tests. In order to facilitate identification of bird species (which can sometime 
vary in quite minor plumage details, especially at certain times of the year), photo¬ 
graphs should be taken according to the guidance given in Sect. 3.2. In the event of 
positive results, further examination of such photos can provide additional informa¬ 
tion on the age and sex of the bird, in addition to proving the identity of the species 
beyond doubt and thus allowing the case to be correctly put into context. To facili¬ 
tate this, each individual bird should be given a code that is used on the cloacal and 
oropharyngeal swabs taken, and this code should be on a piece of card that is visible 
in each photograph taken. 

Especially related to sampling in the vicinity of outbreaks, it is desirable to col¬ 
lect a range of contextual information so as to better understand the viral epidemiol¬ 
ogy of H5N1 HPAI in wild bird populations. Such information should include: 

• date of sampling, clear locational and descriptive data about the catching site, 
ideally GPS coordinates, including habitat description (e.g., lake, river, village 
pond, fish farm, etc.) and distance to human settlement, agricultural land, and 
poultry farms; it may also be useful to include details about the season and relate 
this information to the natural behaviour/cycle of the affected birds, e.g., moult¬ 
ing, premigration, during migration, etc.; 

• record of the numbers of each species of other live birds in the sampling area that 
were not sampled; 

• if available, records of bird movements (arrivals/departures) that occurred at the 
sampling site prior to the sampling; 

• assessment of the numbers of each species of live bird in the sampling area that 
were not sampled but that were showing signs of ill health; and 

• given that birds of some species (such as Mallards-Anas platyrhynchos) can 
occur either as wild birds that are able to move between sites or occur in a feral 
state, habituated to foods provided by humans, distinguishing between these cat¬ 
egories would be useful. Sometimes the presence of unusual plumage patterns, 
indicating domestication, is useful in this respect. 
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Taking photographs of dead birds for identification purposes 
Different bird species are identified by differing characteristics, so it is difficult 
to provide universal guidance applicable in all situations. However, the following is 
a minimum standard that should be followed. All wild birds collected for analysis 
for HPAI should have digital photographs taken as soon as possible after collection. 
The bird should fully fill the photograph and wherever possible include a ruler or 
other scale measure. Each photograph should be taken at the highest resolution pos¬ 
sible and if the camera has a ‘date stamp’ feature then this should be enabled so that 
the image is saved with a time reference - this may help verify the sequence of 
images taken at a site on a day. Images should be downloaded to a computer as soon 
as possible and information about location and date added to the file properties. 
Photographs should be taken of: 

• The whole bird, dorsal side, with one wing stretched out and tail spread and visible; 

• The head in profile clearly showing the beak; 

• Close-up photos of the tips of wing feathers, as these can often determine whether 
the bird is an adult or a juvenile (bird in its first year); and 

• Ideally photographs of both dorsal and ventral views of the bird, as photos of the 
upper and under surfaces of the wing and spread tail will facilitate aging and sex- 
ing of birds 

Any ventral photographs should show the legs and feet (since leg color is often 
an important species diagnostic). If any rings (metal or plastic) are present on the 
legs, these should be photographed in situ as well as recording ring details. Any 
conspicuous markings/patterns should also be photographed. 

At certain times of the year, such as late summer (July-August in the northern 
hemisphere) many waterbirds, especially ducks and geese, undergo moult and can 
be especially difficult to identify by non-specialists. At such times clear photo¬ 
graphs are especially important to aid identification of (duck) carcasses. The patch 
of colour on the open wing (called the “speculum”) is often especially useful. The 
identification of young gulls at any time of the year is also difficult and typically 
they will also need to be photographed and identified by specialists. 

Photographs should be retained, linked to an individual specimen, at least until 
laboratory tests are returned as negative for avian influenza. A unique code or refer¬ 
ence number that is the same as the code or reference number of any samples taken 
from the birds should be visible in each photograph so as to link samples and pho¬ 
tographs. Photographs can be used immediately if identification of the species of 
bird is in any doubt, and for subsequent checking of the identification if necessary. 
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